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The  one-dinensional  fastytine  averaged  Hamiltonian  is  derived 
in  a  free  electron  laser  (FEL)  for  electrons  passing  through  a 

constant  viggler  and  a  radiation  field.  The  exact  unperturbed  orbits 

/ 

vithout  sidebands  ^stfe  obtained  for  all  particles  and  arbitrary 
separatrix  width  8y  Integration  in  action-angle  variables  of  the 
linearized  Vlasov  equation  vith  perturbing  sidebands  over  the 
unperturbed  orbits  yields  the  sideband  gain  for  both  trapped  and 
untrapped  particles.  The  unperturbed  distribution  fg  is  in  an 
adiabatic  equilibrium  vith  the  nain  signal  field.  It  is  found  that 
upper  and  lover  sidebands  that  are  symmetric  relative  to  the  FEL 
frequency  have  opposite  growth  rates.  There  is  no  differentiation  in 
the  magnitude  of  the  gain  betveen  upper  and  lover  sidebands.  The 
stability  is  determine^  b^^the  sign  of  dfQ/d<^y  i.e.t  the  relative 
population  of  oscillation  quanta  «  bounce  frequency  around 

resonance.  The  shear  dco^/dJ,  vhere  J  is  the  action  variable,  is 
stabilizing  and  distributions  vith  gradients  dfQ/dJ  localized  near  the 
separatrix  have  the  minimum  growth  rates.  The  structure  and  scaling  of 
the  unstable  spectrum  are  different  from  previous  results  obtained  for 
electrons  localized  at  the  bottom  of  the  ponderomotive  veil. 
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I.  INTRODUCTION  AND  SUMMARY 


Th«  growth  of  parasitic  aodas  at  frequencies  near  the  aain 

signal  frequency  during  high  pover  PEL  operation  was  theoretically 
1  2 

predicted  *  in  early  1980's.  Since  then  there  has  been  aaple 
3  4  5  6 

numerical’  and  experimental *  evidence  of  sideband  excitation  in 
constant  viggler  FELs.  Unstable  nodes  in  variable  viggler  FELs  have 
also  been  observed  in  simulations^  and  recently  in  experiment1*^. 
Sidebands  degrade  the  sain  signal  efficiency  and  optical  quality  by 
channeling  a  considerable  fraction  of  the  pover  into  parasitic 
frequencies.  The  perfornance  of  the  airrors  in  an  oscillator  can  be 
haraed  fron  the  Modulation  of  the  wave  envelope  caused  by  the 
sidebands.  Last,  but  not  least,  interaction  aaong  nearby  sidebands 
above  a  certain  anplitude  nay  lead  to  chaotic  particle  notion,  loss  of 
trapping  and  incoherent  radiation. 

The  above  have  stinulated  a  considerable  anount  of  theoretical 
work  focused  on  sideband  growth.  Sinple  one-dinensional  configurations 
that  are  analytically  tractable  have  been  used  to  nodel  the  situation. 
Two  lines  of  approach  have  been  considered.  The  single  particle 
picture  regards  the  particle  trajectories  as  functions  of  the  initial 
conditions  and  conputes  the  gain  by  ensemble  averaging  over  initial 
distributions  .  The  alternative  approach  assumes  some  adiabatic 
librium  between  the  particles  and  the  main  signal  and  examines  the 
stability  of  the  perturbations  induced  by  the  sidebands,  solving  the 
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kinetic  equation  *  .  Because  of  the  equilibrium  assumption  the 

kinetic  eethod  is  nore  appropriate  for  PEL  operation  as  an  aaplifier. 

In  both  treatnents  so  far,  analytic  results  have  been  obtained 
only  for  particles  localized  near  the  bo t toe  of  the  ponderoaotive 
veil.  This  inplies  the  following  linitations:  The  sideband  spectrua 
becomes  discrete 


•r  ±  (kr/k^nm^O), 


Vs,* 2v 


where  *#^(0)  is  the  bounce  frequency  at  the  bottom  of  the  ponderoaotive 

veil,  kf  ky  are  the  radiation  and  viggler  vave  numbers  respectively 
2  2  -1/2 

and  Y  -(1-v  /c  )  .  The  contribution  from  untrapped  particles  and 

Z  Z 

trapped  particles  avay  from  the  bottom  is  neglected.  The  effect  of  the 
shear  dw^/dJ,  where  the  action  J  parametrizes  the  distance  from  the 
centre  of  the  separatrix  J*0,  is  lost.  Finally  no  predictions  can  be 
made  about  the  saturation  level  of  potentially  unstable  nodes. 

Here  canonical  formalism  is  introduced  by  expressing  the 
unperturbed  particle  orbits  in  terms  of  action-angle  variables.  The 
unperturbed  orbits  are  of  course  the  fast  time  averaged  "synchrotron'' 
oscillations  of  the  electrons  in  the  potential  veil  formed  by  the 
combined  action  of  the  viggler  and  the  radiation  signal.  The  perturbed 
kinetic  equation  is  solved  in  action  space,  starting  from  an 
equilibrium  extending  over  all  trapped  and  untrapped  electrons.  It  is 


found  that:  0 .  \  — 

(a)  the  spectrum  becomes  continuous  replacing  *^(0)  by  m^(J) 
in  Eq.  (1).  The  modes  located  at  the  peaks  of  the  unstable  spectrum 

1  icutlon _ _ 

grow  faster,  emerging  as  the  discrete  spectrum  that  is  observed  in - — 
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simulations.  5^1  button/ 
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(b)  Nor*  than  one  group  of  particles  are  in  resonance  with  a 
given  sideband  frequency-  a  through  different  haraonics  of  their 

9 

bounce  frequency  and  contribute  to  the  grovth  rate. 

(c)  The  shear  da^/dJ  is  stabilizing.  Distributions  vith 
gradients  dfQ/dJ  localized  near  the  separatrix  are  found  to  have  the 
aininua  grovth  rates  because  of  the  high  shear  there.  This  type  of 
distribution  is  relevant  to  PEL'S  vith  tapered  vigglers. 

(d)  The  gain  is  proportional  to  [df(J)/d«^(J)),  the  relative 
population  in  oscillation  quanta  around  resonance,  in  agreeaent  vith 
the  quantua  aechanical  interpretation. 

(e)  Upper  and  lover  sidebands  located  syaaetrically  around  the 
aaln  signal  frequency  have  opposite  gains  (coaplenentary  stability). 
Therefore  one  node  is  alvays  unstable.  There  is  no  stable  distribution 
fQ(J)  except  the  trivial  one  dfQ/dJ*0. 

(f)  For  any  saooth  distribution,  of  finite  dfQ/dJ,  electrons 
at  the  bottoa  of  the  veil  have  a  negligible  effect  on  stability. 

(g)  Previous  results,  finding  lover  sidebands  having  an 
inherently  larger  gain  than  upper  sidebands,  are  relevant  only  to  the 
Halting  case  of  a  singular  6- function  distribution  f0(J)*&(J).  This 
case  is  unrealistic  because  a  vide,  saooth  initial  distribution  in 
action  fQ(J)  corresponds  to  even  an  ideal  cold  beaa  distribution  in 
aoaentua  fQ(p)*S(p-p0). 

The  nonlinear  saturation  levels  for  the  unstable  nodes  and  the 
aaplitude  for  stochastic  transition  can  also  be  derived  froa  our 
foraalisa  and  vill  be  addressed  in  future  vork. 
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II.  COMPUTATION  OP  THE  GAIN 


Ve  consider  s  linearly  polarized  viggler  of  constant 
wavelength  2n/kw  and  constant  aaplitude  Ay  and  plain  nonochroaatic 
waves  for  the  sain  signal  and  the  sideband.  The  total  vector  potential 
is  given  by  A«Ay(z,t)ey  with 

Ay(s,t)  -  Awcoskwz  ♦  Arcos(krz-*»rt)  +  Agcos(kgz-ngt)  (2) 

where  the  subscripts  v,  r  and  s  stand  for  viggler,  aain  radiation  and 

sideband  respectively,  A  ,  A  ,  are  constant  and  A >>A  .  Our 

*  s  r  3 

Haailtonlan  is  one-diaensional  non  autonoaous  because  the  canonical 
aoaenta 


Px  *  px’  Py  "  Py*  l  V2*1)  (3) 

are  constants  of  notion  that  can  be  conveniently  set  to  zero.  Ve 
switch  to  the  tine  averaged  coordinates  t  •<?  >  and  i-<z>  through  the 
appropriate  canonical  transfornation,  choosing  a  generating  function 
that  kills  the  fast  oscillating  terns  due  to  the  wiggling  notion  (see 
Ref.  13).  Introducing  the  nev  variables 

p  »  pz  /  (hr+  hy),  ♦  -  <kr+  ky)z  -  *»rt 

and  keeping  only  linear  terns  in  the  snail  anplitude  a  ,  the  averaged 

8 

interaction  Haailtonlan  becones 


(A) 


H(P,t|t)-{M2+(kr+kv)2P2-«¥«rco*#-«wagco8(f-4#t)}1/2 

H«r«  tine  is  normalized  to  nr“2,  length  to  kr“*,  velocity  to  c,  mass 
to  ng,  aj-eA^/a^c2  and  M2-l+(ay2+ar2)/2.  Although  kr«#r«l  in  the  above 
units  ve  will  vrite  them  explicitly  to  avoid  confusion.  The  tera  Sg  in 
Bq.  (A)  signifies  the  Doppler  shifted  frequency  departure  of  the 
sideband  froa  the  aain  signal 

«,  -  <ky  /kr)  <*g-  ar).  (5) 

The  terns  proportional  to  a  a  and  a  a  are  the  ponderoaotive 

w  r  w  s 

potentials  caused  by  the  beating  together  of  the  viggling  notion  vith 
the  signal  and  the  sideband  field  respectively. 

When  a#-0  the  Haailtonian  HQ(P,f)  is  exactly  integrable.  The 
unperturbed  trajectories,  shovn  in  Pig.  1(a),  are  given  by 

H0(P,P)  -  R 

vith  the  constant  R  deterained  froa  the  initial  conditions 
R*HQ(P0,40).  These  trajectories  take  the  simplest  possible  fora 
expressed  in  terns  of  the  action  angle  variables  (J,9) 

J(t)  -  J  -  const. , 
o 

e(t)  -  e0  ♦  eebCJ)t,  (6) 

n^J)  -  dHQ(J)/dJ, 

where  J  is  the  area  in  phase  space  enclosed  by  an  orbit  divided  by  2 it. 
In  the  above  variables  the  vidth  of  the  separatrix  v  is  given  by 


¥mj  vith  Jf  the  value  of  action  at  the  separatrix.  The  relation 

aep  aep 

betveen  net  (Jt9)  and  old  (Ptf)  variables  is  expressed  by 


,  ^fj1/2^  (tj(X)-(l-X2)E1(X)l, 

J  '{  *fcF2x . 


sin 


.  X  snI^1(X)ej, 

^  “l  sn(^B^(l/X)0] , 


X2<1 


X2>1 

X2<1 

X2>1 


(7) 


vith  Ej.Bj  the  coaplete  elliptic  integrals  of  the  first  and  second 

kind  and  an  the  Jacobi  elliptic  sine  function.  The  trapping  paraaeter 
2 

X  is  given  by 


X2-  (F+G)/2G, 

r  -  (kr+kv)2(K2-M2),  (8) 

G  -  avar[(kr+kv)2-  1J, 

2  2 

and  ve  have  X  <1,  X  >1  for  trapped  and  untrapped  particles 

2 

respectively.  The  three  constants  of  the  notion  X  ,  R  and  J  are 
autually  dependent  and  any  one  of  then  defines  uniquely  a  trajectory. 
The  bounce  frequency  is  found  froa  Eqs.  (6),  (7)  and  (8) 
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(9) 


V0)ajro’ 

V0>ijW’ 

with 

V0) .  .  !<vm»w. 


X2<1 


X2>1 


Ve  now  pcrforn  a  linear  expansion  of  the  Haailtonian  Eq.  (4) 

in  the  -snail  sideband  aaplitude  a  followed  by  a  decoaposition  of  the 

perturbing  tera  cos[*(9)-&  t]  into  harnonics  of  the  synchrotron 

s 

oscillation  to  obtain 


H(J»0)  -  Ho(J)  +nJj{  h*(J)cos(n©f«st)  +  h*(J)cos(n©-«st)  >, 
hn<J>  '  JjW  F.  <frJ>-  <10) 


The  coefficients  QR  (J)  are  obtained  by  contour  integration  in  the 
coaplex  plane  using  the  doub 
elliptic  functions. 


«£•- 


<*Dn— p 


n 

? 


BJ(X)  l-(-q) 


'n 


n»2X2  n 


lfd/X) 


(  l-oZn 


periodicity  properties  of  the  Jacobi 

JiBj(XK 

0 

q  -  exp 

X2<1 

.  1 

x2<i 

1  .  2n 

1+q 

■  J,  q  ■  exp  ^ 

Ejd/x);’  x  <l 

either  0 

or  "  while  the  aaxiaua  occurs 

2 

at  the  separatrix  X  -1.  Qn(J)  becoaes  progressively  saaller  with 


increasing  n  everywhere  except  in  the  vicinity  of  the  separatricies. 
That  liaits  the  nuaber  of  haraonics  to  be  considered  for  instabilities 
localized  away  froa  the  separatrix.  The  perturbed  equations  of  notion 
in  the  standard  fora  read 


3t  J  -nfonChn<J>sln<n®f*st>  +  -*gt)}, 

dh+  dh~ 

af® -  VJ)  +n!o{ajncos<ne+*st)  ♦  ajnc08<ne-4st))* 

The  gain  gg  for  a  given  sideband  is  deternined  by  the  energy 
balance  equation 


d  2 

a?  v 


2*sV 


lr4nic 


<jy«s>  +  CC)* 


(12) 


The  bracket  <  >  denotes  the  average  of  a  quantity  <f>«l/T  f(t)dt  over 
the  fast  tine  scale  T  -  2,l/kwve  K<  Tb  *  2n/<4b’  and  iy  ls  the 
transverse  current 


jy(*.t)  -  enfa  JJvy4f ( J , 0; t )dJd0, 


(13) 


with 


vy  -  r'lpy  -  -(e/rac)Ay(z,t). 

For  snail  sideband  signal  &f  is  the  pertrurbed  distribution  f^ 


under  the  Haailtonian  flow  Eqs.  (11) 


(14) 


“i 
5T  + 


8fi  dj  «o 

VJ)55 - dt  33“’ 


where  fQ  is  the  equilibrium  distribution  along  the  unperturbed  orbits 
dfo/dt*0.  Any  distribution  fQ(J)  depending  on  J  alone  is  invariant 
under  Bqs.  (6) 


dJdio 
Jf  dJ  * 


0 


thus  fQ  is  not  uniquely  defined. 

Solving  Eq.  (14)  for  fj,  then  substituting  Eqs.  (10)  and  (13) 
into  (12)  and  time-averaging  one  obtains  the  final  result 


g* 


-*4 


2  2 

Wt  11  “nb 

*s  8rt“t 


E 


l°n<J  >  I'  («.]  (i\\1 

|Ho(Jn)*Tr|WJ  'j  WJ  Jj  ' 

n  n 


(16) 


vith  the  bean  plasma  frequency.  The  superscript  ♦  or  -  corresponds 
to  upper  tts>«»r  end  lover  «»a<«r  sideband  respectively.  JR  is  given 
implicitly  by  the  resonant  condition 


*nW  *  <kv/kr>  <VV’  <17) 

as  illustrated  in  Fig.  1(b).  The  sum  over  n  on  the  right-hand  side  of 
Eq.  (16)  includes  the  contribution  from  all  resonant  groups  of 
particles.  The  action  Jn  in  Eq.  (17)  labels  the  orbit  having  the  n-th 
harmonic  of  the  local  bounce  frequency  in  resonance  vith  the  sideband. 
The  gain  is  determined  by  the  slope  of  the  distribution  function  fQ 
near  these  resonant  orbits  J«J  . 


III.  RESULTS  AND  CONCLUSIONS 


The  opposite  signs  in  the  right-hand  side  of  Eq.  (16)  denote 

that  upper  and  lover  sidebands  located  symmetrically  around  the  aain 

signal  |« +-*_|  -  |»-m  |  have  opposite  growth  rates,  g+/g"  «  -1  + 
s  r  s  r 

2 

0( &*/(•_) .  The  physical  reason  is  connected  vith  quant ua  aechanical 

considerations.  The  practical  iaplication  is  that  one  sideband  is 

alvays  unstable  except  in  case  of  trivial  equilibrium  dfQ/dJ-0.  It 

also  reveals  that,  once  destabilized,  upper  and  lover  nodes  have 

grovth  rates  of  equal  magnitude  and  are  equally  dangerous.  Ve  clarify 

that,  in  general,  the  opposite  signs  do  not  imply  that  all  upper 

sidebands  have  the  sane  kind  of  stability,  opposite  to  the  stability 

of  all  lover  sidebands.  Depending  upon  fQ  stability  may  change  sign 

betveen  tvo  upper  (or  lover)  frequencies  because  the  slopes  dfQ/dJn 

change  as  the  location  of  the  resonant  JR's  shift  vith  «g. 

For  a  nonotonic  distribution  fQ  ve  observe  that  trapped  and 

untrapped  particles  yield  opposite  contributions  to  a  given  mode 

because  dct^/dJ  changes  sign  across  the  separatrix.  If  trapped 

particles  are  stabilizing  untrapped  are  destabilizing  and  vice  versa. 

High  shear  is  stabilizing  tending  to  reduce  the  magnitude  of  the  gain 

g  .  This  is  expected  as  the  number  of  resonant  orbits  is  inversely 
8 

proportional  to  |dtt^/dJ|.  Shear  tends  to  infinity  near  the  separatrix 
thus  the  nodes  generated  by  electrons  orbiting  near  the  separatrix, 
corresponding  to  small  |«g-«r|  «  (kr./kw)n«jj(0),  have  the  smallest 

grovth  rates.  Nodes  coming  from  electrons  near  the  center  Jn«0  vith 


11 


|  *  (kr/kw)n«|j(0)  have  also  snail  grovth  rates  as  the 

coefficients  Qn(J)  tend  to  zero  there.  In  particular,  given  any  snooth 
distribution  fQ(J),  electrons  exactly  at  the  botton  of  the  veil  have  a 
null  contribution  to  the  instability.  It  takes  a  singular  distribution 
of  infinite  gradient  at  J-0  to  create  an  instability  from  electrons 
exactly  at  the  center  of  the  separatrix.  The  shear  is  generally  higher 
for  untrapped  particles  thus  the  trapped  ones  usually  dominate  the 
instability. 

Although  Eq.  (16)  vas  obtained  on  purely  classical  arguments 
it  nevertheless  admits  the  correct  quantum  mechanical  interpretation. 
Expressing  (df0/dJ)(do^/dJ)~*  as  d^/dw^  ve  observe  that  the  gain  gg 
is  proportional  to  the  difference  in  population  of  oscillation  quanta 
between  the  energy  levels  across  the  resonance.  A  more  detailed 
explanation  for  the  opposite  signs  in  Eq.  (16)  is  given  in  Ref.  13. 

The  normalized  gain  g  /«  is  plotted  against  the  percentage 

S  l 

mismatch  (u  -»  )/«  for  both  upper  and  lover  sidebands  in  Figs.  2  and 
s  r  r 

4.  The  contribution  up  to  the  third  harmonic  n<4  in  Eq.  (16)  is 

included  in  these  plots.  The  parameters  chosen  correspond  to  a  viggler 

-4 

wavelength  X  -3  cm,  a  -5,  main  signal  strength  a  -5x10  beam  energy 

W  W  S 

of  11.43MeV  (0-0.999,  y-22.37)  and  current  density  j«100  A/cm2  (beam 
10  -3 

density  6.25x10  cm  ).  The  equlibrium  distribution  is  a  Gaussian 
?  1/7  7  2 

fo(J)«(l/2nIr)A/*exp(-(J-J0)  /2D  j  centered  halfway  inside  the  island, 

J0«v/2,  and  of  width  D  equal  to  half  the  separatrix  distance  D-v/2, 

vith  v-J  .  Ve  plot  the  contribution  of  only  the  fundamental,  n-1  in 
sep 

Eq.  (16),  in  Fig.  2(a)  adding  the  first  harmonic  n«2  in  (b)  the  second 
harmonic  in  (c)  and  the  third  harmonic  in  (d).  New  unstable  bands 
emerge  vith  each  harmonic,  while  the  gain  for  already  unstable  bands 


is  Modified.  For  example  ve  observe  tvo  upper  and  tvo  lover  unstable 
bands  in  (b)  but  only  one  upper  and  three  lover  bands  in  (d).  Ve  find 
the  contributions  from  higher  than  the  third  haraonics  n>4  generally 
negligible.  The  upper  frequency  c#  >«  and  the  lover  frequency  «  <« 

parts  of  the  unstable  spectrua  come  from  the  regimes  of  negative  and 

positive  slope  df^da^  respectively,  shown  in  Fig.  3.  The  lover 
sideband  grovth  is  peaking  at  frequencies  corresponding  to  no>h(Jaav) 
vith  J||ax  the  value  maximizing  df  /da^fJ).  The  peaks  for  the  upper 

sideband  grovth,  hovever,  do  not  occur  at  J-0  that  minimizes 

dfo/du^(J),  but  at  J  halfvay  inside  the  negative  slope  regime.  This  is 
because  Q  (J)  and  consequently  g  are  zero  at  J«0,  shoving  the 

o  s 

negligible  contribution  from  electrons  at  the  bottom  of  the  veil.  The 
analogous  effect  in  plasma  physics  is  the  elimination  of  the  thermal 
effects  vhen  pL  goes  to  zero.  In  any  case  the  most  unstable  modes  are 
far  from  the  fequencies  lws-wr  |"(k  /k  )n«L  (0)  pointed  by  the  arrovs 
in  Fig.  2(d). 

Since  sidebands  can  not  be  shut  out  completely  it  remains 

debatable  vhether  a  distribution  function  can  be  tailored 

experimentally  to  minimize  their  grovth  rate.  From  the  previous 

discussion  a  flat  distribution  inside  the  trapped  regime  vith  sharp 

gradients  localized  at  the  separatrix  seems  the  appropriate  choice. 

Instabilities  vill  then  localize  near  the  separatrix  and  the  gain  vill 

be  suppressed  by  the  strong  shear.  To  check  this  ve  plot  the  grovth 

rates  from  Eq.  (16)  for  tvo  types  of  distributions  fQ< J) s  (i)  Tvo 

2  1/2  2  2 

Gaussians  fQ(J)-(l/2iiD  )  exp(-J  /2D  )  centered  at  the  centre  of  the 
island  and  of  characteristic  lengths  D  equal  to  half  the  island  vidth 
D-v/2  in  Fig.  4(a)  and  one  island  vidth  D*v  in  Fig.  4(b).  (ii)  Tvo 


"step-like"  distributions  of  the  fore  f0(J)»(l/ed>)*/*exp(-(J/ed))N) 
vith  N-16.  Selecting  o^(N/N-l)^/N  places  the  sharp  gradient  at  J-D  and 
ve  plot  the  case  D»v/2  in  Fig.  4(c)  and  D*v  in  Fig.  4(d). 

Coe paring  Fig.  4(a)  to  Fig.  4(b)  and  Fig.  4(c)  to  Fig.  4(d)  it 

is  seen  that  the  growth  rates  betveen  similar  types  of  distributions 

tend  to  decrease  as  the  location  of  the  eaxieue  gradient  (df  /dJ)  v 

o  ux 

approaches  the  separatrix.  In  both  cases  there  is  nore  than  one  order 
of  magnitude  reduction  in  the  gain  by  shifting  the  maximum  gradient 
position  from  D-w/2  to  D«v.  Because  fQ  was  chosen  monotonic  in  all 
above  plots  and  because  it  was  limited  to  trapped  particles*  dfQ/du^ 
preserves  sign  and  only  lover  nodes  are  unstable.  The  spectral  width 
of  the  unstable  regimes  is  reduced  vith  a  parallel  increase  in  the 
maximum  gain  as  one  goes  from  the  Gaussian  type  to  the  step-like  type 
of  distributions.  Also  the  distance  of  the  sideband  frequencies  from 
the  main  signal  decreases  by  shifting  the  gradient  position  D  closer 
to  the  separatrix.  Distributions  vith  sharp  gradients  at  the 
separatrix  such  as  those  in  Figs.  4(b)  and  4(d)  are  perhaps  more 
relevant  to  the  case  of  variable  viggler  PEL,  vhere  the  "bucket"  of 
the  trapped  particles  is  decelerating  in  phase  space  leaving  the 
untrapped  particles  behind,  and  only  small  diffusion  occurs  across  the 
separatrix  alloving  sharp  gradients. 

The  limit  of  a  8-function  distribution  fo«8(J-JQ),  examined 
13 

elsevhere  ,  is  the  most  unstable  case  but  of  the  least  practical 
interest,  because  even  the  case  of  a  monoenergetic  beam  distribution 
Pz  «po  is  described  in  J-space  by  a  smooth  distribution  fQ(J)  of 
finite  vidth  AJ  (see  Fig.  1). 

Our  calculations  done  for  the  case  of  a  fixed  wavelength 


viggler  can  be  easily  extended  to  the  case  of  a  variable  viggler  PEL 
provided  that  the  saae  adiabatic  assuaptions  hold.  Only  the  functional 
relations  Eqs.  (7)  between  the  action-angle  variables  (J,0)  and  the 
coordinates  (P,t)  need  to  change.  The  derivation  of  the  growth  rate, 
performed  in  J-space,  is  indepedent  of  the  trans format ion  J(P,t)  and 
the  result,  Bq.  (16),  stands  as  is.  It  is  the  high  gain  of  tapered 
wiggler  FBLs  that  seems  to  challenge  the  adiabatic  approach  to  the 
problem.  In  this  case  the  change  of  the  signal  amplitude  af(t)  in  time 
together  with  the  depedence  of  the  pseudo-equilibrium  fQ  on  both  time 
and  9  should  be  included  for  a  more  realistic  treatment. 
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Figure  1.  Tine  averaged  notion  without  the  aideheads.  (a)Plots  in 

phase  space  of  the  unperturbed  orbits  H0(P,+)-K.  The 

intersections  with  the  horisontal  line  P>const.  nark  the 

initial  conditions  for  each  orbit,  (b)  The  noraalised  bounce 

frequency  and  the  first  two  haraoaics  as  functions  of  the 

2 

trapping  paraneter  X  (J).  The  intersections  with  the 
horisontal  line  «Mnf-nr)k¥/kr  determine  the  position  JQ  of 
the  resonant  orbits  for  a  given 

Figure  2.  Gain  for  a  Gaussian  distribution  f0(J)«Cexp(-(J-Jo)*/2D*) 
centered  halfway  inside  the  trapped  particle  island  J0«w/2 
and  of  width  D  equal  to  w  where  The  noraalised  gain 

gs/«r  is  plotted  versus  S*/«r  for  (a)  the  fundamental 
contribution  n-1  in  Iq.  (16),  (b)  including  the  first 
harmonic  n-2,  (c)  two  harmonics  and  (d)  three  harmonics.  The 
nost  unstable  nodes  do  not  correspond  to  harmonics  of  the 
bounce  frequency  at  the  bottom  e^(0)  indicated  by  the 
arrows. 

Figure  3.  Plot  of  dfo(J)/de^(J)«df0/dJ(dU|)/dJ)'1  as  a  function  of 
K(J)-K(0)  for  the  distribution  fQ  of  Figure  2.  The  slope 
goes  to  sero  near  the  separatrix  K(JMp)  because  of  the 
infinite  shear  du^/dJ. 

Figure  4.  Normalised  gain  for  aonotonic  distributions  centered  at  the 
bottom  of  the  well  J»0  including  the  first  three  harmonics 
n£t  in  Eq.  (16).  (a)  Gaussian  distribution  of  width  D  equal 
to  half  the  island  width  v,  D»v/2  (b)  Gaussian  distribution 
with  D>w.  (c)  Step-like  distribution  with  D-w/2  and  (d)Step- 
like  distribution  with  D-v. 


17 


XIO 


Flf.  5 


o  1 

<0 

I 


Fig.  4 


